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ABSTRACT 

We address the question of whether or not the properties of all low-luminosity 
flat spectrum radio sources, not just the obvious BL Lac objects, are consistent with 
them being the relativistically beamed counterparts of the low radio luminosity ra- 
dio galaxies (the Fanaroff- Riley type 1 - FR I). We have accumulated data on a 
well-defined sample of low redshift, core-dominated, radio sources all of which have 
one-sided core-jet structures seen with VLBI, just like most BL Lac objects. We first 
compare the emission-line luminosities of the sample of core-dominated radio sources 
with a matched sample of FR I radio galaxies. The emission lines in the core-dominated 
objects are on average significantly more luminous than those in the comparison sam- 
ple, inconsistent with the simplest unified models in which there is no orientation de- 
pendence of the line emission. We then compare the properties of our core-dominated 
sample with those of a sample of radio-emitting UGC galaxies selected without bias 
to core strength. The core-dominated objects fit well on the UGC correlation between 
line emission and radio core strength found by Verdoes Kleijn et al. (2002). The re- 
sults are not consistent with all the objects participating in a simple unified model in 
which the observed line emission is orientation independent, though they could fit a 
single, unified model provided that some FR I radio galaxies have emission line regions 
which become more visible when viewed along the jet axis. However, they are equally 
consistent with a scenario in which, for the majority of objects, beaming has minimal 
effect on the observed core luminosities of a large fraction of the FR I population 
and that intrinsically stronger cores simply give rise to stronger emission lines. We 
conclude that FR I unification is much more complex than usually portrayed, and 
models combining beaming with an intrinsic relationship between core and emission 
line strengths need to be explored. 

Key words: galaxies: galaxies: active - galaxies: radio - galaxies:jets - galaxies: emission 
lines. 



The utility of unified schemes of active galaxies is their sim- 
plicity and their potential for making testable predictions. 
In the case of powerful FR I radio sources some of the initial 
simplicity has been lost because modifications have had to 
be introduced in order to match the predictions to the obser- 
vations (see Urry & Padovani, 1995 for a review). In the case 
of low power radio sources, however, the simple idea that 
FR I radio galaxies looked at down their jet axes become 



1 INTRODUCTION 



BL Lac objects has remained the consensus view (Padovani 
& Urry, 2001). This is perhaps somewhat surprising since 
only about a third of low-luminosity core-dominated radio 
sources (supposedly the beamed counterparts of FR Is) are 
found to be conventional BL Lac objects; most of the re- 
maining two-thirds of such objects have optical classifica- 
tions such as Seyfert-like objects or passive elliptical galax- 
ies (Marcha, et al., 1996). What is the relationship between 
these other types of core-dominated objects and FR Is? 



* e-mail:mmarcha@oal. ul.pt 



We have been exploring the idea that the synchrotron 
cores of all low luminosity core-dominated radio sources are 
intrinsically similar, irrespective of their emission line prop- 
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erties, and that all might be appropriately labeled blazars. 
Such a scheme is consistent with their continuum SEDs 
(Caccianiga & Marcha, 2004; Anton et al., 2004; Anton & 
Browne, 2005). 

In this paper we investigate the emission-line properties 
of low luminosity radio galaxies, both core-dominated and 
lobe-dominated. We first compare the emission-line proper- 
ties of a group of core-dominated objects with those of a 
matched sample of FR I radio galaxies. According to the 
simplest unified scheme, in which there is no orientation de- 
pendence of the line emission, the emission line properties of 
the two samples should be statistically indistinguishable. A 
rigorous exploration is possible because we now have avail- 
able a sample of core-dominated radio sources - the 200 mjy 
sample (Marcha et al., 1996) - for which we have emis- 
sion line luminosities, extended radio luminosities and high- 
resolution VLBI observations. Knowing the extended radio 
luminosities is potentially important because the emission- 
line properties of AGN correlate with radio power (Rawlings 
& Saunders, 1991; Verdoes Kleijn et al., 2002) and it is there- 
fore necessary to choose a comparison sample from objects 
having the same range of intrinsic radio powers. From a uni- 
fied model point of view, the weak extended radio emission 
of the core-dominated objects (i.e. excluding the core and 
one-sided jet emission) should be a good measure of the in- 
trinsic radio power since the symmetry of the emission, and 
its low surface brightness, eliminates the possibility that its 
observed strength is influenced by relativistic beaming. The 
high-resolution VLBI observations are also very important 
since we cannot simply rely on spectral selection to pick 
core-dominated sources since not all radio sources with flat 
spectra have radio structures dominated by cores. Some are 
found to be Compact Symmetric Objects (CSOs) in which 
the flat spectrum does not arise from compact self-absorbed 
core emission. VLBI observations enable us to select with 
confidence the genuine one-sided core-jet objects and dis- 
card the rest. 

A second line of investigation can be adopted. If the 
non-BL Lac core-dominated objects are not relativistically 
beamed, then we are measuring the intrinsic core strengths 
and thus objects with strong cores might be expected to 
produce more line emission. Thus, looking at correlation be- 
tween line emission and radio core strength might be useful. 
Verdoes Kleijn et al. (2002) have done this for a sample 
of UGC galaxies selected without reference to radio core 
strength and find quite a strong correlation. From the tight- 
ness of the correlation they conclude that the bulk Lorentz 
factors lie in the range 2-5 for continuous jets, or <2 for 
jets consisting of discrete blobs. Most of the objects in their 
sample are not core-dominated. It would, therefore, be in- 
teresting to see if other objects with more dominant radio 
cores (in which the effects of beaming would be much more 
apparent if there existed bulk relativistic motions) lie on a 
continuation of their correlation. We can use the data on the 
200 mjy sample and combine them with the UGC results 
for this purpose. 

The structure of the paper is as follows. In Section 2 we 
describe the selection of the core-dominated sample and of 
a comparison sample matched in flux density and redshift 
to the first. In Section 3, we briefly describe some radio ob- 
servations of members of the core-dominated sample that 
were made in order to determine their extended radio flux 



densities. In Section 4 we compare the emission line distri- 
butions of the two samples and also investigate where the 
200 mjy objects lie on the Verdoes Kleijn et al. (2002) cor- 
relation. Finally, in Section 5, we discuss the results, both in 
terms of unified models and in terms of "dis-unified" mod- 
els. Throughout we adopt a value for the Hubble constant 
of 75 km s" 1 Mpc" 1 . 



2 THE SELECTION OF THE SAMPLES 

We start with the 200 mjy sample (Marcha et al., 1996) 
updated by Anton et al. (2004). This sample consists of 
flat-spectrum sources (Sex v~ a with a <0.5) stronger than 
200 mjy at 5 GHz, as recorded in GB6 (Gregory et al., 1996), 
with red optical magnitude brighter than 17. The magnitude 
selection means that out to a redshift of ~0.1 we have an 
essentially volume-limited complete sample. The sample is 
further restricted to declinations > +20° and Galactic lat- 
itudes >12°. Measurements of the_ luminosities of the H a 
lines are available for most objectaD. VLBI maps have been 
made of virtually all the objects at 5 GHz and/or 1.6 GHz 
(Bondi et al., 2001; Bondi et al., 2004; Bondi & Polatidis, 
personal communication) and we have only included in our 
discussion those objects with a core and one-sided jet visi- 
ble in the VLBI maps - future maps may well reveal more 
core-jet objects. We have also restricted our sample for the 
current paper to objects with redshifts <0.1 since in this red- 
shift range we are confident that the sample is complete. The 
redshift restriction also facilitates finding a good compari- 
son sample (See below and Table 1). Spectroscopically, the 
sample has been classified according to the following types 
(also identified in Table 1): PEG - stands for Passive Ellipti- 
cal Galaxies and it refers to sources with weak emission lines 
and a strong galaxy component; Syl,2 - stands for sources 
with strong emission lines in their spectra (1 if broad emis- 
sion lines are seen, 2 if only narrow lines are present); hyb - 
stands for sources with broad but significantly weaker emis- 
sion lines than in the Sy case; BL Lac - stands for sources 
with weak or absent emission lines in their spectra. Seven 
BL Lac objects with continuous optical spectra and no mea- 
sured redshifts from the 200 mjy sample are excluded from 
the sample to be discussed. This should not bias the statis- 
tics because it is almost certain that all these objects are at 
redshifts >0.1, otherwise the host galaxy would be visible 
and its redshift known. 

We use as a basis of our comparison sample the ra- 
dio galaxies found in Abell clusters and studied by Owen 
et al. (1995) and by Ledlow & Owen (1995). These consist 
of objects with 20 cm flux densities >10 mjy and redshifts 
<0.09. H Q luminosities (or limits on luminosities) are given 
by Ledlow & Owen for most of the objects. These sources 
are mostly ellipticals from the spectroscopic point view, al- 
though a small number (less than 13%) has strong emission 
line spectra related to starformation (Owen et al. 1996). 
In our analysis we will use a sub-sample of the Ledlow & 
Owen sample consisting of objects with H a measurements 
and picked to match the 200 mjy objects both in redshift 

I When we refer to H a we implicitly include emission from the 
[Nil] lines with which H Q is usually blended. 
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and in flux densityHJ For each 200 mjy object we have looked 
for a "twin" matched to within a factor of two in 20 cm flux 
density to the 200mJy object's extended flux density (see 
next section) and v2 in redshift (in most cases the matches 
are much closer) , regardless of the strength of the line emis- 
sion. Matches were possible in all but one case. The proce- 
dure used to produce a comparison sample eliminates any 
worry about the effect of the correlation between line lumi- 
nosity and the intrinsic radio power. It also ensures that the 
linear scales of the spectrograph slit, when projected on the 
galaxy, are comparable. The twin for each 200 mjy source 
is given in Table 1. In this process we are implicitly assum- 
ing that the Ledlow & Owen objects are representative of 
the population of unbeamed counterparts of the 200 mjy 
objects. 



3 VLA OBSERVATIONS AND DATA 
REDUCTION 

The best estimate we have of the intrinsic radio power of 
a core-dominated object, where the core and jet emission 
may be Doppler boosted, is the extended radio emission. 
Estimates of the extended emission can be made using the 
FIRST survey (Becker et al., 1995) data but these are avail- 
able for only a fraction of the 200 mjy objects. For this 
reason we have made our own VLA B-configuration obser- 
vations at 20 cm of those 200 mjy objects which do not 
lie in the region covered by the FIRST survey. Each source 
was observed at two widely separated hour angles for a to- 
tal observing time of ~10 min. The observations were cal- 
ibrated and mapped using standard AIPS tasks. The data 
were of good quality and in most of the maps, extended ra- 
dio structure is detectable. The peak flux densities and total 
flux densities were measured with the task IMEAN and the 
extended flux density taken as the difference between the 
two quantities. The results are listed in Table 1. For those 
sources in the FIRST survey, FITS images were retrieved 
and processed in the same manner as our own maps. 



4 THE EMISSION LINE PROPERTIES OF 
CORE-DOMINATED AND 
LOBE-DOMINATED OBJECTS 

4.1 The 200 mJy sample and comparison sample 

We first ask the simple question: are the emission line prop- 
erties of the 200mJy sample sources statistically distinguish- 
able from those of the sources in the comparison sample? We 
choose to compare H a luminosities since these are widely 
available and do not worry at this stage whether we are 
dealing with narrow-line or broad-line emission. In the sim- 
plest version of the FR I unified scheme in which there are 
no hidden emission line regions one would expect the distri- 
bution of emission line luminosities in the 200 mjy sample 

t Since we have good radio spectra for the 200 mjy sources we 
do the matching in extended flux density at 20 cm. The choice of 
wavelength to do the matching does not affect the validity of the 
comparison. 



and in the comparison sample to be indistinguishable (see 
Figure |l|) . 

An advantage of having matched pairs of objects is that 
it gives us a simple way to quantify the probability that 
the distribution of line luminosities for the two samples is 
indistinguishable. In 16 pairs (excluding 3 pairs where both 
members have limits on luminosities and one where they 
were identical) it is found that in 13 cases the member drawn 
from the 200 mjy sample has the more luminous H a line. 
We use the binomial distribution to work out the probability 
that in 13 or more cases the 200 mjy member should have 
the larger line luminosity, on the assumption that the line 
luminosity distributions of both groups are the same. We 
conclude that there is less than 2% probability that the line 
luminosity distribution in the two samples is drawn from the 
same population. 

The confidence in this result could be questioned due 
to two factors: firstly there may be some degeneracy in the 
choice of twins (e.g. more than one possible twin for each of 
sources of the 200 mjy sample), and secondly, there is the 
statistically difficult issue of comparing distributions where 
there is a significant number of upper limits. We decide to 
investigate the consequences of both of these aspects by pro- 
ceeding in the following way: 

• For each 200 mjy sample source we take all the twins 
that fall in the region within a factor of two in 20cm flux 
density and a factor of v2 in redshift; 

• For the 16 sources with 4 twins or more we take the 
line luminosities (the detections and upper limits) and find 
the Kaplan-Meier estimator for the distribution function of a 
randomly censored data by using survival analysis (ASURV; 
LaValley, Isobe and Feigelson, 1990). For the remaining 
sources, we preferentially took the twin with a detection, 
or in case there were only upper limits, the highest of these 
as the line luminosity for the twin. 

• The line luminosity of the twins found in the way de- 
scribed above is then compared to the line luminosity of the 
200 mjy sources. We find that in 17 cases the line luminosity 
of the 200 mjy sources is higher than that of the twin, in 1 
case it is equal, and in 2 cases it is smaller. If we then use 
the binomial distribution to estimate the probability that 
in 17 out of the possible 19 pairs, the line luminosity of the 
200 mjy source is larger than that of the comparison sample 
we find that this probability is <0.0001. 

• Finally, we note that 17/20 sources of the twin sample 
have L(H a ) < 10 33 W, whilst in the 200 mjy sample, only 
3/20 have such low values of line luminosity. 

We therefore think that the difference between the 
strength of the line luminosity of 200 mjy sample and its 
twin should be taken seriously. There is at least one selec- 
tion effect which is likely to mask the true difference be- 
tween core- and lobe-dominated objects when we use our 
two samples. The two samples are cross-contaminated; the 
200mJy sample contains some not very core-dominated ob- 
jects while the Ledlow and Owen sample contains some ob- 
jects which are just as core-dominated as the 200mJy ob- 
jects. For example 0055+300 (NGC315) a giant FR J radio 
galaxy (Bridle et al., 1976) is in the 200 mjy sampleBI. Such 

§ It is expected that some radio galaxies with strong radio cores 
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Figure 1. Distributions of the line luminosity for the 200 mjy sample (left panel), and its twin sample (right panel and marked SI). 
The histograms marked with <, correspond to upper limits. 



a cross-contamination is likely to make it more difficult for 
us to recognize true differences in emission-line properties 
rather than leading to spurious differences. Hence, the de- 
tected difference is likely to be a lower limit on the true 
difference between the line luminosities of the two samples. 

There is one caveat related to the fact that the Led- 
low and Owen objects are selected to lie in Abell clusters of 
galaxies whereas the 200 mjy sources are selected without 
any reference to their cluster environments. It could be that 
the hosts of radio sources in clusters have different emission- 
line properties to galaxies in the field. Such an effect has been 
claimed to exist by Guthrie (1981). We are, however, reas- 
sured that this may not be too important an effect since the 
200 mjy objects appear to be consistent with the correlation 
for non-core-dominated objects found by Verdoes Kleijn et 
al. (2002) between emission line luminosity and radio core 
luminosity (see below and Figure ^). 

4.2 The 200 mJy sample and UGC galaxies 

The result above suggests a simple interpretation which is 
that emission line strength is related to core strength. Such 
a view is supported by the observation of a tight correlation 
of the core emission line strength with radio core strength 
in the UGC FR I radio-galaxies (Verdoes Kleijn et al., 2002; 
Xu et al., 2000; Figure 2). Given that most of the residual 
dispersion in the correlation could be due to measurement 
errors and/or to radio core variability, that leaves little room 
for the effects of beaming on the core strengths. Variability 
is relevant for narrow-line emission because the Narrow Line 
Region (NLR) probably has a size of ~tens to hundreds of 
parsecs and thus its observed strength reflects the core activ- 
ity integrated over times-scales of tens, perhaps hundreds, of 
years. Hence, one spot measurement of radio core strength 
may not be a good indicator of the integrated core activity. 



will find their way into the 200 mjy sample simply because of 
random errors on the 5 GHz and 1.4 GHz flux densities used to 
define the sample 
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Figure 2. Distribution of the line luminosity vs. the radio core 
luminosity for the UGC FR Is of Verdoes Kleijn et al. (filled 
squares) the 200 mjy sample (open symbols; triangles are the BL 
Lacs). 



To investigate if the correlation holds for objects con- 
tained in the 200 mjy sample, we include these sources 
with the UGC objects of Verdoes Kleijn et al. in a plot 
of the VLBI core luminosity (L rc ) against the line lumi- 
nosity L(Ha) (Figure where we separate the BL Lacs in 
the 200 mjy from the remaining sources by plotting them 
as open triangles). Because the 200 mjy sources have more 
dominant cores than the UGC objects, they might be ex- 
pected to display the effects of beaming more strongly. (We 
note that the comparison sample is not included in this plot 
since they lack VLBI measurements). This should manifest 
itself in two ways: (i) as an increase in the dispersion in 
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the correlation, and (ii) in the more core-dominated objects 
falling oil average below the line established for the UGC 
objects. 01 It is clear that the 200 mjy objects join relatively 
smoothly on to the existing UGC correlation. 

The correlation coefficient for the UGC+200 mjy ob- 
jects (excluding the 6 BL Lacs in the 200 mjy where 4 have 
upper limits to the line luminosity) is 0.72 corresponding 
to a probability of there being no correlation of <Cl%. The 
linear least squares fit gives a slope of 0.89±0.1. There is an 
increase in the scatter going from the UGC to the 200 mjy 
objects but no obvious sign in the decrease in slope which 
would be expected if the radio core emission is beamed 
and the line emission isotropic. In fact, the BL Lacs of the 
200 mjy sample seem to fit within the general relationship 
without significantly increasing the dispersion of the dis- 
tribution. This is confirmed by the generalised Spearman's 
correlation test of 0.76, when we use survival analysis to 
include the upper limits of the line luminosity of the BL 
Lacs. This value means that the probability of no correla- 
tion between the two quantities is less than 1%. The result- 
ing linear regression when all the sources are considered is 
L[H a ) oc (0.76±0.09)L rc . 

We emphasise that the correlation between radio core 
and line luminosity is not induced by redshift. In fact, the 
partial correlation coefficient between the two quantities, ex- 
cluding the effect of redshift is 0.62 for the UGC and 200 mjy 
sources (where the six BL Lacs were not considered) which 
gives a probability << 1% to the hypothesis that there is 
no correlation between the two luminosities. 



5 DISCUSSION 

The basis of our discussion is that the more core-dominated 
sources have significantly stronger line emission than the 
lobe-dominated sources. This shows up in the comparison 
of the 200 mjy sample objects with the matched sample 
objects and in the strong correlation between radio core lu- 
minosity and emission line luminosity. What, if anything, 
does this tells us about FR I unification? We will start by 
looking at things from two opposite viewpoints, one unified 
and one "dis- unified" . We then discuss if there is a middle 
way. 

5.1 The unified interpretation 

If we adopt a strict "orientation is everything" point of view, 
we conclude that in some FR Is, optical emission line regions 
must often be hidden from the observer, and (statistically) 
become more visible when looking close to the jet axis. We 
note that in most cases, but not all, the line emission we 
see in the core-dominated objects of the 200 mjy sample 
is not nuclear broad-line emission (which is believed to be 
present in most FR lis but sometimes hidden from view 
by a molecular torus), but most often narrow-line emission 

> We note that the plotted line luminosities are slightly differ- 
ent in the two samples. In the case of the UGC objects the line 
emission refers to the nuclear H a + [NII] emission measured with 
HST, whereas in the case of the 200 mjy objects we are using the 
H a + [NII] flux integrated over a wider spatial extent. 



which originates at radii more than a few parsecs from the 
nucleus of the galaxy. More specifically, there are four ob- 
jects, 0125+487, 0321+340, 1646+499 and 2116+81, in the 
200mJy sample in Table 1 which are classified as having 
broad emission lines and all four have high line luminosi- 
ties. Although better (higher resolution) spectra would be 
required to provide narrow line luminosities, visual inspec- 
tion of the available spectra shows that all the broad line 
sources have substantial narrow line components which are 
probably sufficiently strong by themselves for the galaxies 
to be called Seyfert 2s. This means that in nearly all cases 
it is the narrow line emission, which is thought to originate 
more than a few parsecs from the nucleus, that would have 
to be hidden. This type of narrow-line emission would not 
be affected by a torus but might be partially hidden by, for 
example, extinction in a kind of dust and gas disks seen 
m many FR I radio galaxies (de Ruiter et al., 2002)0. Ad- 
ditional evidence for off-nuclear gas and dust is presented 
by Quillen et al. (2003) who argued, on the basis of 3CRR 
radio-galaxy Spectral Energy Distributions (SEDs) , that the 
nuclear regions of some FR I radio galaxies suffer from opti- 
cal extinction. Further evidence comes from recent results by 
Wills et al. (2004) on 13 low luminosity sources of the 2 Jy 
sample. In particular, these authors found that the average 
[OIII] line luminosity in BL Lacs is significantly larger than 
that of FR Is, thus supporting the view of extinction in the 
latter type of sources. Since a factor of a few in orientation- 
dependent extinction in the inner kpc would be enough to 
account for the excess of line emission in the core-dominated 
objects with respect to those in the comparison sample, we 
think this is just about a viable scenario. 

5.2 A dis- unified interpretation 

An alternative interpretation is to say that the observations 
rule out FR I unification or, at least, force us to exclude 
from any such schemes the core-dominated synchrotron jet 
sources which are not recognizably BL Lacs. This would not 
be a very tidy state of affairs, but it could be the truth! Such 
a view, one might argue, is the most logical interpretation of 
the tight correlation of the core emission line strength with 
radio core strength in the UGC FR I radio-galaxies (Verdoes 
Kleijn et al., 2002), and the fact that the 200 mjy objects lie 
on the same correlation (Figure 2). To explore this correla- 
tion further we have added to the objects plotted in Figure 
2 those FR1 radio sources from the 3CR studied in Cao & 
Rawlings (2004) - see Figure ^ where we have excluded the 
3CR FR I galaxies that were in common with those already 
present in the UGC sample. We emphasize that no attempt 
has been made to match the 200 mjy sample objects in 
extended luminosity. It is clear that the core luminosities 
cover a similar range to those seen in the 200 mjy sample 
and the line luminosities are consistent with objects from 
this additional sample participating in the same correlation 
as the objects from the other two samples. Simple beaming 
schemes, even those incorporating orientation dependent ex- 
tinction of emission line regions would not predict this. 

II We note that Chiaberge et al. (2002) argue that the combina- 
tion of optical and UV HST data, and radio core emission for 
FRIs is not consistent with the standard thick torus framework. 
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Figure 3. Distribution of the line luminosity vs. the radio core 
luminosity for the UGC FR Is of Verdoes Kleijn et al. (filled 
squares) the 200 mjy sample (open symbols; triangles are the 
BL Lacs), and the 3CR FR 1 galaxies of Cao & Rawlings (filled 
triangles). The radio core flux densities were converted to 1.4 GHz 
using a radio spectral index a r = 



One might argue that because the slope of the 
core/emission line correlation is less than unity, a viable 
hypothesis is that there is an underlying linear correlation 
between H a luminosity and intrinsic core strength, and all 
the scatter in Figure 2 could be attributable to the effects 
of beaming. In other words, beaming should leave H a un- 
changed but move the core luminosity towards the right in 
Figure 2. We have tested this idea in two ways. First we look 
for a correlation between the magnitude of the deviation of 
the core luminosity from the underlying linear (slope unity) 
correlation and the core dominance parameter R (L rc /L ex t 
at 1.4 GHz). Both should be orientation indicators in the 
beaming scenario. We find absolutely no correlation (see Fig- 
ure ^) , which reinforces the view that the evidence for radio 
core beaming effects amongst these objects is very weak. The 
second approach is to adopt the traditional unified scheme 
view (e.g. Orr & Browne, 1983) that there is a tight relation 
between the intrinsic core strength and that of the extended 
lobe emission. In this case, the measured core-dominance (R) 
compared to what would be seen if the object were viewed 
in the plane of the sky (R c ) is a direct measure of the core 
Doppler boosting. Based on models which use population 
statistics to constrain beaming model parameters, current 
estimates for R c are about 0.01 (Wall and Jackson, 1997; 
Jamrozy, 2004). We use these estimates, and the observed 
R values, to de-boost each object. The resulting plot of line 
luminosity against de-boosted core luminosity is shown in 
Figure ^ where symbols are as before. The original correla- 
tion is destroyed. 

One additional thing we have tried, still within the con- 
text of there being a linear correlation between line lumi- 
nosity and intrinsic core strength, is to estimate how much 
orientation dependent extinction would be required to just 



Figure 4. Distribution of the residuals of the line luminosity - 
radio core luminosity correlation (see text for details) against the 
radio core parameter R for the UGC FR Is of Verdoes Kleijn et 
al. (filled points) and for the 200 mjy sample (empty symbols; 
triangles are BL Lacs) 

counterbalance the effect of beaming and restore Figure 5 
to look like Figure 3. The result of our exercise, however, 
is that it would require one to move most of the 200 mjy 
objects in Figure 5 by around two orders of magnitude in 
line luminosity which is much larger than current estimates 
of the extinction suffered by the cores of most FR I radio 
galaxies would allow. For instance, based on the study of the 
UV emission of the nuclei of 3CR galaxies, Chiaberge et al. 
(2002) deduce a median value of Av = 1.3 for the FR Is in 
their sample. 

5.3 Reconciling the interpretations 

We are faced with a dilemma. The evidence for relativistic 
bulk motion in the cores of BL Lacs is extremely strong 
and hence the core emission we observe in these must be 
highly orientation dependent because of Doppler boosting. 
There is also direct evidence from superluminal motion that 
there are highly relativistic jets at the pc-scale in FR1 radio 
galaxies (Giovannini et al., 2001). The case for unification 
between BL Lacs and some, or all, FRls is thus very strong. 
On the other hand, we have shown that BL Lacs are not 
the only low-luminosity core-dominated radio sources with 
synchrotron jets (Marcha et al, 1996; Bondi et al., 2001; 
Anton et al., 2004, Bondi et al., 2004), and is very tempting 
to hypothesise that these non-BL Lacs have relativistic jets 
too and that they fit into a single unified picture. We note 
for instance that two PEGs (1241+735 and 2320+203) of the 
200 mjy sample are completely indistinguishable from BL 
Lacs not only in pc-scale morphology, but also in the level of 
polarization detected in the pc-scale jet (Bondi et al., 2004). 
However, in this paper we find that, despite the similarity in 
radio properties, the non-BL Lacs have on average stronger 
emission lines than the FR1 radio galaxies. Furthermore, 



© 2004 RAS, MNRAS 000, 1-?? 



Radio galaxy unification 7 



36 



34 - 



□ A 

"a T □ 



32 - 



20 22 
log I. (W/Hz) (debeamed) 



21 



Figure 5. Distribution of the line luminosity vs. the radio core 
luminosity for the UGC FR Is of Verdoes Kleijn et al. (filled 
points) and for the 200 mjy sample (open symbols) after the core 
luminosity has been 'de-boosted' (see text for details). 



the strong correlation between observed core strength and 
emission line strength, both for the UGC sample alone and 
the UGC sample plus 200 mjy, is prima facie evidence for 
a tight relationship between the intrinsic core strength and 
the emission lines that are somehow excited by this core 
activity. Is it possible to fit all we know about these objects 
into a single framework? 

The simplest interpretation is to postulate an orienta- 
tion dependence of the emission lines which makes them 
appear systematically stronger as the angle to the line of 
sight decreases. This could go some way to account for why 
stronger emission lines go with stronger cores but fails on 
three counts: 

(i) It does not explain why the BL Lacs, which are be- 
lieved to be viewed at the smallest angles, do not have strong 
lines. 

(ii) It does not account for the very tight correlation be- 
tween radio core strength and line emission in the UGC sam- 
ple which has been selected on extended optical and radio 
properties, and which should not introduce any orientation 
bias. 

(iii) The amount of extinction (~ 5 m ) required to restore 
the correlation in Figure [| seems much too large when com- 
pared with the recent measurements by Chiaberge (2002) 
for FR I radio galaxies. 

Any complete unified model must have an intrinsic cor- 
relation between core strength and line emission built in 
from the start. This may be possible. In a subsequent paper 
we will explore a model based on the following ideas: 

• Jets have a highly relativistic spine and a slower sheath 
(Laing & Bridle, 2004). 

• The relativistic spine of the jet is where the characteris- 
tic BL Lac emission comes from and is usually only dominant 
when viewed at small angles to the line of sight. 



• It is the energy dumped into the slow sheath that drives 
the production of emission lines and produces the majority 
of the observed core emission at large viewing angles. 



6 SUMMARY AND CONCLUSIONS 

We summarize our results as follows: 

• We have compared the H a emission properties of the 
200 mjy sample of low-luminosity core-dominated sources 
with those of a matched sample of FR I radio galaxies and 
find that the former have significantly stronger emission 
lines. Either unified models based on beaming have to be 
modified to include some orientation-dependent extinction, 
or the dispersion in intrinsic core strengths is so large that 
beaming is not the dominant factor. 

• From combining observations of UGC galaxies with 
those of our 200 mjy sample we see that the observed core 
and emission line strengths in radio galaxies are strongly 
correlated, with little evidence for the kind of behaviour one 
might expect if the radio emission from the cores of most of 
the core-dominated objects were beamed. The evidence for 
beaming amongst the many BL Lacs may be strong but it 
is far from strong amongst the many objects with virtually 
identical radio properties. 

We conclude that FR I unification is much more com- 
plex than usually portrayed and models combining beaming 
with an intrinsic relationship between core and emission line 
strengths need to be explored. 
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Table 1. The 200mjy sample and the comparison sample. Cols: (1) Name; (2 and 11) Extended flux density at 20cm; (3 and 10) rcdhsift; 
(4) spectroscopic type: Syl,2 stands for Sy-type 1,2; PEG stands for Passive Elliptical Galaxies; hyb stands for hybrid; BL stands for 
BL Lac; (5) log of line luminosity; (6) log of extendend luminosity; (7 and 12) equivalent width of H a + [Nil] in A; (8) SED type taken 
from Anton et al., 2004: SPL for Steep Power-law, bPL for broken power-law and PL+IR concave spectrum with one or more bumps; 
(9) Name of twin (B1950). Flux densities are in mjy, radio luminosities in W Hz^ 1 and line luminosities in W. 
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